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ABSTRACT 


Long-duration surface missions to the Moon and Mars will require habitats for the 
astronauts. The materials chosen for the habitat walls play a direct role in the 
protection against the harsh environments found on the surface. Choosing the 
best materials, their configuration, and the amount required is extremely difficult 
due to the immense size of the design region. Advanced optimization techniques 
are necessary for habitat wall design. Standard optimization techniques are not 
suitable for problems with such large search spaces; therefore, a habitat design 
optimization tool utilizing genetic algorithms has been developed. Genetic 
algorithms use a " survival of the fittest" philosophy, where the most fit 
individuals are more likely to survive and reproduce. This habitat design 
optimization tool is a multi-objective formulation of structural analysis, heat loss, 
radiation protection, and meteoroid protection. This paper presents the research 
and development of this tool. 

Keywords: lunar habitat, design, optimization, analysis 



INTRODUCTION 


The nation’s Vision for Space Exploration calls for a human return to the Moon 
by 2020 and a later potential human mission to Mars. These future missions are 
expected to involve much more impressive activities than those of the Apollo 
program. Some missions may last several months, while others may last up to 
600 days (1). These long duration surface missions require large outposts to 
accommodate living quarters (habitats) for the astronauts as well as indoor 
laboratory facilities. Transporting the materials required to build the necessary 
habitats is expensive (2). The greatest impediment to extended human presence 
on the Moon or Mars is the threat posed by the harsh environments found on their 
surfaces and in transit. The lunar environment is much different than the 
terrestrial environment; many of the load conditions on the Moon are extremely 
severe (2). 

Crew habitats require impact resistance, thermal insulation, radiation protection, 
and structural support. The materials chosen for the habitat walls play a direct 
role in protection against each of the mentioned hazards. The cost of launching 
payload from Earth is approximately $10,000 per pound. Through the use of in- 
situ resources and terrestrial materials a habitat can be designed to protect against 
the hazards on the surface of the Moon and Mars. 

Choosing the best materials, their configuration, and the amount required is 
extremely difficult due to the immense size of the design region. Clearly some 
sort of optimization method would be extremely beneficial. 

OPTIMIZATION 

There are several optimization techniques that could be used, each having a 
number of advantages and disadvantages. Standard optimization techniques such 
as hill climbing, Newton’s method, etc., are not suitable for problems with the 
large search space that will be encountered when designing a surface habitat, 
because they seek only local optima (3). 

Genetic algorithms are search algorithms based on the principles of evolution and 
natural genetics. They use a “survival of the fittest” philosophy, where the fittest 
individual is most likely to survive and reproduce, therefore transmitting its genes 
to the next generation (4). Genetic operators - crossover and mutation - are used 
to create the succeeding generation from the members (strings) of the current 
population. Each individual is assigned a fitness value based on how well it 



satisfies the objective function, which is generated through engineering analysis 
of the candidate design. The probability of a design being propagated to the next 
generation as well as the frequency of occurrence of the design in the next 
generation is directly proportional to the fitness of the individual (3). 

BACKGROUND 

Objective 

The purpose of this work is to develop a cross-disciplinary surface habitat 
optimization and analysis tool that can determine the best combination of chosen 
materials for the walls of a surface habitat. The tool optimizes the habitat 
materials for minimum upmass including the analysis of thermal losses and gains, 
structural integrity, meteoroid impact, radiation shielding, and material upmass. 

A pictorial depicting the various objectives is shown in Figure 1. 

Meteoroid Impact 

A meteoroid is a “naturally occurring solid body, traveling through space, which 
is too small to be called an asteroid or a comet” (5). It is estimated that 
meteoroids of approximately milligram mass strike lunar structures yearly. 
Velocities of meteoroids at the moon range from 13-18 km/s. Meteoroids of 
about IQ' 6 g can generate craters of 500pm in metals. The impacts of larger 
meteoroids are rarer; however, they pose a more significant threat. A meteoroid 
with a mass of about l g can produce a crater of centimeter scale in metals (5). 

Thermal 

The thermal stability of a lunar habitat is of extreme importance. A lunar habitat 
experiences drastic thermal cycling loads during a lunar day. Temperature ranges 
can vary 200K on the surface of the moon; therefore, it is necessary to design a 
structure that can sustain such dramatic thermal cycling loads. Active heating and 
cooling is required for lunar habitats and a reduction of such a need will reduce 
the upmass of the required equipment. Therefore, thermal losses are calculated to 
determine the suitability of a given design. 

Structural 

Structural integrity of the lunar habitat is necessary to ensure that the habitat will 
remain intact. A lunar habitat must be able to withstand the residual loads from 
the habitat itself and the pressure differential of the “shirt sleeve” environment 
and the vacuum of space. All of these loads are considered when evaluating the 
structural integrity of the lunar habitat. A safety factor of 1.1 of the yield strength 
was used to analyze a loaded habitat. 



Radiation 


The radiation environment found in interplanetary space consists primarily of 
energetic ions and the galactic cosmic rays (GCR). Further sources of radiation 
are solar particle events (SPE). These are infrequent events that occur with large 
coronal mass ejections (CME) and result in a large amount of particles, mostly 
protons, moving through the solar system (6). 

The presented habitat design tool calculates the dose equivalent on the inside of 
the habitat using simple material density and composition correlations. 

Materials Database 

The materials database contains all the necessary material properties for each of 
the analysis modules. Currently there are 24 materials in the database including 
aluminum, polyethylene, Mylar, lunar regolith cement, polyetheretherketone 
(PEEK), stainless steel, Kevlar and glass composites, rubber, and ethylene-vinyl 
alcohol copolymer (EVOH). 

CODE DEVELOPMENT 

The code was written primarily in MATLAB. Each analysis module was written 
as a separate subroutine. The optimization code is a genetic algorithm written by 
the authors at NASA Marshall Space Flight Center. 

Meteoroid Impact 

The meteoroid impact analysis is centered on the Fish-Summers single plate 
penetration equation (7): 


Ti 0.352 0.167 /jr o \ 0 875 \ 

P = Km p p p (Fpcos P) (1) 

where P = depth of penetration (cm), K = material constant, m p = mass of the 
projectile, p p = density of the projectile, V p = velocity of the projectile, and (3 = 
angle of impact. 

The projectile density and velocity were assumed to be 0.5 g/cm and 20 km/s, 
respectively. The velocity is conservative, and the density is an estimate. A zero 
degree angle of impact was used to be conservative. 

The material constant can be approximated for each material with the following 
equation (8): 
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where s = percent elongation and p T = target density. 


The projectile mass was determined (9) to be 0.22g by iterating between Equation 
3 and 4 (10): 


F = c 0 {(c,w 0306 +c 2 ) 438 +c 2 (m + c 4 m 2 +c 5 w 4 ) 0,36 +c 6 (m + c 1 m 2 ) 0,85 }(3) 

where F = number of particles/m 2 /year of mass, m, or greater, and the constants 
are listed in Table I. 


P* = e~ F * A *‘ (4) 

where PNI = probability of no impact, F = number of particles/m 2 /year of mass, 
m, or greater, A = area (m 2 ), and t = time (years). 

To be conservative, the area was calculated based on a 32m outer radius, 15 year 
time, and 95% probability of no impact. 

Only the outside layer of the habitat wall was analyzed by the impact code. This 
was chosen, because averaging all of the material constants (K) for each layer of 
the multi-layered wall would have resulted in less accurate penetration depths 
than by only analyzing the outer layer. Furthermore, this forces a tough outer 
layer to protect the habitat, which is most realistic. 

Thermal 

Steady State Numerical Thermal Model 

A significant amount of research has been done to develop an analytical model of 
the multi-layer insulation heat leak of the habitat. Analysis of the habitat begins 
by the application of the first law of thermodynamics to the control volume. 

This is a passive system, which means that no additional heat comes in or is taken 
away by thermal control system. A depiction of this system is shown in Figure 2, 
and the corresponding energy balance equation is the following: 

K +E g + E out = E sl 


(5) 



where E m is the energy into the system from outside sources, E out is the energy 
loss from the system. E g is the thermal energy generation and E st is the rate of 
energy stored within the system. 

The core problem of cooling or heating lunar habitats is to combat the solar 
heating by direct and indirect means: directly from the sun and indirectly from 
lunar reflection. The effective radiation heat transfer equation for the sun on the 
habitat is the following: 


QradSun ~ ° ‘ 4 Hab ' ( 6 ) 

where a is absorbtivity, A Hab is incident surface area of the habitat, and G s is the 
solar constant. 

The effective radiation from the habitat takes the form of: 


QradHab 
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where e is the emmissivity, a is the Stephen-Boltzman constant, A is the surface 
area of the habitat, and T s is the surface temperature of the habitat. 

Calculation of the surface temperature of the habitat requires the determination of 
the roots of a forth degree polynomial. Equating the conduction and radiation 
heat transfer rates and the polynomial is shown in Equation 8. 
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Thermal Model 

A finite element thermal model of the thermal behavior of the lunar habitat was 
created in ANSYS, and an example result is shown in Figure 3. The model uses 
steady state analysis for a 2D cross section of the habitat. 



Structural 


A structural model was created and written in ANSYS. This structural model 
includes all loading conditions including pressure differential, thermal loading, 
and material weight in a 1/6 g environment. An example ANSYS plot is shown in 
Figure 4. 


Radiation 

The radiation analysis is a rule-of-thumb, level-zero analysis. It converts the 
thickness of each layer into a corresponding thickness of aluminum and then 
calculates the dose equivalent in the middle of the habitat. 

The 1989 Solar Maximum radiation environment is used in this analysis, because 
it is assumed to be a near- worst-case solar particle event. The radiation model 
only accounts for protons, but a 5 rem dose is added to account for any secondary 
neutrons created from the primary proton bombardment. 

Materials Database 

The materials database is contained as a matrix that is called when the tool is 
started. It contains all the necessary material properties required for each of the 
analysis modules. 

Optimization 

A genetic algorithm (GA) was written specifically for the lunar habitat wall 
optimization problem. This GA uses a combination of crossover, mutation, and 
attrition to search the design region for the lunar habitat optimization problem. 

The design region is defined as the area containing all of the possible solutions to 
the given problem. A design region of this stature is difficult to search. For 
example a 3 -layer habitat wall created from 20 possible materials all with varying 
thickness would have the following number of possible solutions: 

(number of layers)* (possible material types) (p0SSlble th,ckness) = 320°°. 

However, a GA can successfully find a “good” solution within this range. 

OBJECTIVE FUNCTION FORMULATION 

The habitat optimization tool described in this report attempts to optimally 
combine several aspects pertinent to a winning lunar habitat design: provide 
meteoroid protection, minimize thermal losses, maintain structural integrity, and 



provide radiation protection. Each of the mentioned measured parameters is 
independently calculated based on the material and geometric framework of the 
candidate habitat design. The difficulty lies in the weighting of the individual 
objectives in question. Should they all be weighted equally and how is equality 
measured? These questions must be addressed to folly optimize a realistic lunar 
habitat design. 

The weightings of the objective functions are simply a subjective measure the 
designers ordered priority. It may be the choice of the designer to emphasize the 
minimization of the active thermal control required to maintain the “shirt sleeve” 
environment for the astronauts or the meteoroid and radiation protection. 
Regardless of the designers choice it is necessary to scale the objectives such that 
the weightings are initially equal. Recall the various objectives are measured in 
heat loss, deflection, penetration depth, and radiation dose. To equally weight 
these parameters a series of 60 analyses were performed, the objective constraint 
values were changed sequentially to equalize the differing objective parameters. 
The final values are listed in Table II. 

These values are simply an equal weighting representation of the individual 
objectives. At this point the designer has the option of emphasizing a single 
objective over another. Examples of designer weighted objectives are given in 
the next section. 


EXAMPLE 

An example optimized habitat configuration is shown in Figure 5. This 
configuration was selected with the weightings listed in Table III. The result is a 
habitat created primarily from lunar cement. This is expected since there is such a 
high emphasis on upmass. 


CONCLUSIONS 

A cross-disciplinary surface habitat optimization and analysis tool has been 
developed that can determine the best combination of materials for the walls of a 
surface habitat. The tool optimizes the habitat materials for minimum upmass 
including the analysis of thermal losses and gains, structural integrity, meteoroid 
impact, and radiation shielding. The tool presented in this paper provides a 
preliminary optimized habitat to give a design team a starting point for more 
detailed engineering and analysis. The focus of the optimization can be changed 
by varying the weightings of each of the analysis constraints. This adds a great 



deal of flexibility for the user and allows the design to be focused for each 
application. 

There is an evident value in seeking optimal habitat wall configurations. An 
optimal wall will provide mass savings, radiation protection, impact protection, 
insure a reliable base for future missions, and provide a safe haven for astronauts 
against the harsh environment of the lunar surface. 
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TABLES 


Table I. Constants for Equation 3. 


Constant 

Value 

co 

3,1 56x1 0 7 

Cl 

2.2x103 

C2 

15 

C3 

1.3x10'* 

C 4 

10“ 

C5 

To 27 

C6 

1.3xl0' 16 

C7 

“To 5 


Table II: Final equa 

ly scaled constraint factors. 

Constraint 

Equation 

Constraint Factor 

Upmass 

upmass*upmass_factor 

1 

Heat Loss 

Q * thermalfactor 

11.68 

Structural 

(5 + sf)*structural_factor 

872227.76 


P*penetration_factor 

95233.43 

Radiation Dose 

Rem*radiation_factor 

51064.29 


Table III: Weightinj 

2 ,s for example optimization. 

Constraint 

Weighting 


Upmass 

6 

Heat Loss 

0.3 

Deflection 

0.1 

Penetration Depth 

0.3 

Radiation Dose 

0.3 
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Figure 1: Multi-objective parameters for the lunar habitat design tool. 
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Figure 3 : Example FEA thermal model. 







